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Abstract 

The paper extends the classical result on the convergence of the 
Krawtchouk polynomials to the Hermite polynomials. We provide the 
uniform asymptotic expansion in terms of the Hermite polynomials. 
We explicitly obtain expressions for a few initial terms of this expan¬ 
sion. The research is motivated by the study of ergodic sums of the 
Pascal adic transformation. 
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1 Introduction 


Let 0<p<l,q = l —p and let A be a positive integer. The non-normalized 
Krawtchouk polynomials of variable x can be defined by the identity 


Kn{x,p,N) 


2^1 


—X, —n 
-N 


1 

p 


( 1 ) 


where x and n are in {0,1,... N} and 2F1 is the Gauss hypergeometric 
function. The normalized Krawtchouk polynomials are usually dehned as 
follows: 

kiP\x,N) = (-pP A„(x,p, A). (2) 
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The second argument N is usually omitted, so we write kl^\x) instead of 
{x, N). They form an orthogonal system on the discrete set {0,1, 2,A^} 
with weight function 

ATI N-x 

o{x\ = - 

’ r(i + x)r(iv +1 - x) 

and orthogonality relation 

^kf\x)kf\x)p{x) = { .VpqySij, z,j = 0,l...,iV. 

x=0 ^ 


Krawtchouk polynomials satisfy the following Rodrigues-type formula 
(see [1], Section 2, (22a)) 



(_g)n A^(p(x)x^^) 
n! p{x) 


(3) 


where A/(x) = /(x + 1) — /(x) and y- = y[y — 1)...(?/ — A: -|- 1). Formula 
([3]) is also often taken as the dehnition of the Krawtchouk polynomials, see 
e.g. [6]. Expression ([3]) also shows that the Krawtchouk polynomial kn\-) 
can be considered as an analytic function on [0, A]. 

Finally, it was recently found in |3j that the Krawtchouk polynomial 
{—2p)'^Kn{k,p, N) has a natural interpretation as the ergodic sum along 
tower TN,k of the Pascal adic transformation for the orthogonalize^ Walsh- 
Paley function where the sum of binary digits in the binary representation 
of t G N equals to n (see [3] for details). 

There is a classical result on convergence of the (properly renormalized) 
Krawtchouk polynomials to the Hermite polynomials (see El): 



with X = Np+ {2NpqY/‘^x and Hn{x) = (—The result was 
extended by Sharapudinov in [^, where he obtained the asymptotic formula 


{2NpqTrn\y^^{Npq) "'^^p{x)e^^^‘^klf\x) = 


with X = Np + {2Npqy/‘^x, n = x = 0(n^/^). 

^with respect to the Bernoulli (p, q) measure 
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Let ^4 be a positive real. In this paper we are interested in the higher order 
uniform in u G [—Ay/N, A-\/N\ asymptotic expansion of the Krawtchouk 
polynomials of the form: 

M 

k(P\x) = + ^ (g) 

j=0 

where x = Np + v, n = 0(1), MgNU{ 0} and [t] gives the integer part 
of a real number t. Our main result is the uniform asymptotic expansion of 
p{x)kn\x) in terms of the Hermite polynomials is stated by Theorem[2] Our 
approach is based on the result by V. V. Petrov (j5]) extending the Local 
Limit Theorem (LLT). 

We are especially interested in the expression for the first non-constant 
term Cj{v) in expansion Q, due to the study of ergodic sums of the Pascal 
adic transformation in [3]. It turns out that the value of j depends on the 
parity of index n. In Corollary [1] under the additional assumption v = 
we obtain elegant explicit expressions for c\{v) and C2{v) for odd 
and even values of n correspondingly. 

Recently, there is a considerable interest in the asymptotics of the Kraw¬ 
tchouk polynomials, when the parameter N grows to infinity (see e.g. jT] and 
references therein). In paper [T] authors considered (among many others) the 
case when x = 0(1) while N oo. Nevertheless they ruled out the case n ~ 
Np. We mention, that the self-duality relation Kx(n,p, N) = Kn{x,p, N) 
implies that the one-term asymptotic expansion follows already from (|3]) in 
this case. 

The author is grateful to Profs. A. A. Lodkin and A. M. Vershik for their 
advices and kind support. 

2 Main results 


Using identity A®x— = we rewrite Rodrigues formula ([3]) as follows 


p{x)klPHx) = ^A-[pix)x^] = ^ 


ni 


ni 


k=0 


^ ) A’"-V(x) A^(x + n-k)^ 


{-qT 

n! 


n 


E 



(n — k) 

ky —^A 


n—k 


p{x). 


(7) 


Next we consider the term A^p(x), s > 0, separately. We introduce h- 
step forward difference operator Ahf{x) = f{x + h) — f{x), A))/(x) = 


3 






1 


A/i(A^ ^/(x)),n > 2. Following |6| we let h to be equal to 
can writ^ 

AV(x) = A^/)(x(x)), 


V27V^’ 


thus we 


( 8 ) 


where x = Np + {2Npq)^^‘^x. In fact, asymptotic relation (jj]), as well as 
(j5|), follow^ from the LLT applied to p{x) and the mean-value theorem 
~ h"-f {x + nh9), 6 E (0,1), combined with the proper estimation 
of the residual term (see details in El). In order to obtain higher order 
approximation we consider function p{x) as a probability of x successes in a 
sequence of N independent (p, g)-Bernuolli trials and use Theorem 16 from 
§3 of the work [5] by V. V. Petrov for the Bernoulli distributior0: 


Theorem 1. Let M he a nonnegative integer and cr^ = PQ- The following 
asymptotic expansion holds uniformly in x such that Np+ {2Npq)^/‘^x E Z; 


(1 -h ( y/Np{x) - 


M ~ i 
-2 qv[,x, 




E 

u=Q 


iV!2/2 


= o 


1 


jYM/2 J ’ 


(9) 


where i = Np + (2Npq)^/‘^t and functions q^ are defined as follows: 


1 '' I 

ux )=n — 


Tm+2 


A km'.\{m + 2)\a^+^ 


( 10 ) 


with 'ji, i > 0, being cumulants of the Bernoulli {p,q) distribution, and 
the summation in the right-hand-side taken over all nonnegative solutions 
{ki, ^ 2 , ■ ■ ■, kiy) of the equation ki -\- 2k2 + • • • + I'kiy = v such that s = 
ki k2 ku. 


To obtain the higher order approximation for A^ we use the formal rep¬ 
resentation A/j = — 1 where D = -^ is the difference operator which 

yields 

CXD 

A^ = (e'“°-l)' = (11) 

i=s 

^Note that A is forward difference in x, while Ait is h-step forward difference in x. 

^Other approaches are based on the convergence of a the difference equation, which 
polynomial solutions define the Krawtchouk polynomials, to the differential equation defin¬ 
ing Hermite polynomials, see details in e.g. [4], or on the convergence of generating func¬ 
tions, see [3- 

"^In book [5] a slightly different definition of the Hermite polynomials was used: 
Henf*) = (—l)"e“^ They are related with polynomials Hnix) by the equality 

He„(a:) = 2-"/2Htt(^"). 
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For any nonnegative integer K and any analytic real-valued function / 
we can truncate the series and simply write 


K 

= E + o(h^+*+'). 

i=0 


( 12 ) 


Coefficients could be found by means of the Multinomial theorem as 
follows: 


O-sj 



(13) 


where s and j are nonnegative integers and the summation in the right- 
hand-side of (|13p is taken over all nonnegative solutions (ki,k 2 , ■ ■ ■ ,kj) of 
the equation ki -|- 2 /c 2 + ■ ■ ■ + jkj = j satisfing ki + k 2 + ■ ■ ■ + kj = s. In 
particular we have 



as,2 


s(3s -|- 1) 
24 


and so on ... 


(14) 


Theorem [T] suggests that in order to obtain asymptotic expansion for 
Af^p{x) we need the asymptotic expansions for Qu{x)), v > 0. We 

denote by bu,s coefficients ^^,/ 2 +s) 0 ) arising in the right- 

hand-side of formula (fTOll . Thus we have q^ix) = ^ b^^sHu+ 2 s{x)- We denote 
by gu,rix) the following expression e* qu{x). For Hn{x)) 

we have the identity 




dx 




dx^ ^ dx"- 


e-"') = (-l)^dd„+,(x), 


which implies 


(15) 


~9uA^) = = Y,{-lYb,,,H,+2s+r{x), 

where limits of summation in the right-hand-side are the same as in (|10p . 
Let ^ be a positive real and r be a nonnegative integer. In Appendix we 
show that we can differentiate asymptotic expansion m and the estimate 
of the residual is valid uniformly for anij^ x from the set [—A, A]: 


dx‘ 


:An p{x{x)) = 




ira 


M 

E 

u=0 


9u,r{x) 

iVi^/2 


-h O 


jYM/2 


(16) 


®For s = 0 we use a convention = f{x). 

®Note that Theorem [T] states this for r = 0 and values of x from the discrete set only. 
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It is instructive to compare the above result with Theorem 7, Section 6, 

in [5]. 

Remark 1. Since parameter M in Theorem [1] and formula (|16p is an 
arbitrary positive integer, we can write residual term as Q^ jv(m+i)/ 2 ^ instead 

of ^ • 

Together with expansions (1161) and ([12]) and taking into the account that 
h ~ this results the following asymptotic expansion for A®p(x): 


K 


^s+i 


VNAi p{xix)) = ^ + 0 (iV-(^+^+i)/ 2 ) 


i=0 


M 


’ ^ /)—n j/—n 


i=0 12=0 

/ K M-i _ , ■. 

" ^- - 9u,r\^) 


'' ^ 2=0 l /=0 




+ 0(Ar-(^+^+i)/2) 


where K = M + 1, x = 0(1) and x{x) = Np + {2Npq)^^‘^x. This expression 
is written in fact in terms of the Hermite polynomials, but it also could 
be rewritten using the parabolic cylinder functions Dn{x) due to their close 
relation with the Hermite polynomials: Dn{x) = G 

NU {0}. 

K K-l-i - X 

We denote by il^f{x) the sum ^ ag^i Yh ^‘'/v 02 where coefficients 

i=o ’ v=0 

agj, with j and s G N U {0}, are defined by expression (1131) . 

Using representation ([7|) together with (|17p we obtain our main result 
providing asymptotic expansion of the Krawtchouk polynomials in terms of 
the Hermite polynomials: 


Theorem 2. Let M and n be nonnegative integers and A be a positive real. 
Let X — Np = {2Npqy/'^x. Then the following asymptotic expansion holds 
uniformly in |x| < A: 


p{x)k^P\x) = 


(-qY 


\/2'kN a n\ 


M 

E 

A:=0 


n 


(n —fc) 


n- 


{x+n-k) 2 )• 

(18) 


Proof. Since h ^ N we see that if^{x) = 0{N and 

(x + s)~ = O(N^) and therefore (x + ^)~i^) = 


6 


















Formula (jl8l) follows directly from Rodrigues formula d?]) by changing 
by their approximations ^ (x) obtained in (|17p . 


□ 


In the special case M = 0 Theorem [2] reduces to the well-known result: 



n! 


+ o(iVt), 


stated above in (HI) (we used p(x) = , ^, e -|-o(l) from Theorem [T] here). 

VinNa _ 

In the general case expression (IlSp . in fact, provides expansiorlO ([6]) (i.e. 
its right-hand-side comprises only powers of N and v = Xy/2Npq, not of 
^/N h ^). To see this note, that depending on parity of index n Hermite 
polynomials Hn contain either odd or even powers only, expression m 
does not change parity of indexes of Hermite polynomials while (|15p implies 
expansion (11211 also preserves parity of powers. 

In the special case of M = 2 and v = x — Np = we are going 

to obtain a more natural expression then (slightly cumbersome) expression 
(fTHp suggested by Theorem [2l 

First of all note that for the Hermite polynomials we have the following 
representations (see e.g. 0): 


H2iix) = i-iy2\2l - 1)!!(1 + 

H,i^,{x) = (-l)'2'+i(2/ + l)!!(x + 

where I is a nonnegative integer. Additional assumption x —0 provides 
asymptotic expansions: 

H 2 i{x) = {-iy2\2l - 1)!!(1 - 2lx^) + o{x^), 1 = 1,2..., 

H 2 i+i{x) = {-iy2^+\2l + ly.lx + o{x^), / = 0 , 1, ... 

For any nonnegative integers I, m and C we have the following asymptotic 
expansion for the falling factorial (m + C)- as m —>■ oo: 

(m + Cy- = m} + {iC - (20) 

^Of course, it is sensible to take M <n. 
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Corollary 1. For any sequence e(iV) such that lim £{N) = 0, the following 

N^oo 

asymptotic expansions hold uniformly in v satisfying |?;| < e{N)N^/^ : 

where ti = 6{p — ^)(4/ — 1) and t 2 = {I — 1) (l + 4/ + (16/ — 5)pq) and I G N; 


k^iUNp+ v) = i-l) 


{21 - l)\\{pqN)^ Al{p-\)+2,v 
(2/)! 3 


+ o(iV'), 


where / G N U {0}. 


Proof. We denote by 'ijj{x) the function . Let k he a nonnegative 

integer. Applying Theorem [1] for M = 2 and identity (|15p . we obtain 




73 Hs+k{x) 
23/20-3 3!v/]v 


74 


H^+kix) + ^{^) HQj^k{x) 


+- 


+ O 


1 

iv/’ 


4- 4!iV 

where 73 = pq{l — 2p), 74 = —pq{6pq— 1). Under the additional assumptiorj^ 
V := X — Np = o(Ay3/3) we get a simpler expression: 




V \ ^ ^ _ ^-pq- hv{p - q) J_\ 


( 21 ) 


In the same manner for odd-order derivatives we have 

.s... 36/n^-|-riu-|-r2 




36pqN 



where ri = 6(1 — 2p){2l -|- 3)(2/ -|- 1), T 2 = {21 + l){2l -|- 3)(1 -|- / — (1 -|- Al)pq), 
and for even-order 


j21 


i^{x)^p{x) = {-iy+\2l + l)\l2^ 


V2 , , (l-2p)(2/ + 3), , ,1, 

vwf ^^ 


with I G N. After some algebra using expansions (I19p and (1201) and ignor¬ 
ing terms of the order lower than we obtain from (|18p the required 

asymptotic expansion (and it naturally depends on the parity of n). □ 

®The slightly weaker assumption v = 0{N^^^) would already need additional term 
i {pqyk'N'^ in the right-hand-side of ((21]). 
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Taking higher-order asymptotic expansions in (I19p and (|20p and taking 
appropriately high value of M one can easily obtain other values of Cj, j > 1, 
from expansion Q. 

We conclude this paper with a slight modification of Corollary [1] for the 
function Kn{x,p, Ni), with Ni = N — i, v = x — Np = and i = 0(1). 

Corollary 2. For any sequence e{N) such that lim £{N) = 0, the following 

N^oo 

asymptotic expansions hold uniformly in v satisfying |u| < e{N)N^^^: 


K 2 i{x,p,N-i) 


[ - 1 )!!%v + ipf+ ii{v + ip) + i 2 ^\ ^ 

V p) V j, 


whereii = 6{p—^){4l — l) andi 2 = {l — l){l+4:l + {16l — 5)pq)—9pq{i + 2l — l) 
and I E N. 


K 2 i+i{x,p,N-i) 



{21 + 1){21 - 1)\\ 4l{p 

F+i 


5 ) + 3{v + ip) 
3p 




where / E N U {0}. 

Proof We have Kn{x,p, Ni) = Kn{Np + v,p, Ni) = Kn{Nip + vi,p, Nf) = 
hn\NiP + vi,Ni), where vi =v + ip. Using (l20l) together with 
expression (A^ — i)~^ = {N)~^ -|- il {N)~^~^ + 0{N~^~‘^) we easily obtain the 
required asymptotic expansions from Corollary [TJ □ 

The uniform asymptotic expansions from Corollaries [Hand [2] are obtained 
under the assumption |n| < We can also use the preceding asymp¬ 
totic expansions for |r;| < e{N)N^/‘^, ^{^) 0. However, if = 0{v) we 

have to multiply the estimations of the residual terms of these asymptotic 
expansions by y/N. 


3 Appendix 

In this section we prove formula (|16p . We denote by Sj^{ff) such function 
that 

y/Np{i)=e^Nii)(^^{^)^ 

2 ^ - (S) 

where 4>^ {^) = J e~^ ^ polynomials q^, are defined in Theo- 

rem 1 and, as above, s = Np-\-y/2Npqs. We denote by 4>^(f) the derivative 
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d‘ iM 




M 


(0 = (- 1 ) 


^ • Using the integral Cauchy formula 


we obtain 


i /=0 


(x — 0 ®+^ 


s' 

27ri 


/ 


, f {y/NpiO - 4>^{i))di 

(x — 


(x — 




c. 


(x - 


where C^; is a closed contour around point x. Somewhat arbitrary we set 
to be a unit circle around point x: 


Cx = {^\^ = X + e*^, p G [0, 27r]}. 


Lemma 1. Let M be a nonnegative integer and A be a positive real. Then 
for X G [—A, A] and ^ Cx we have the uniform estimate: 

Proof. The Stirling expansion for the Gamma function (see e.g. book [8], 
p. 34, Example 1.2., or |9], p. 83) could be written as follows: 

Inr(z) = iz-^) ln(z) - z + + F^iz) + Oi\z\-^^-^), (22) 

where | 2 ;| oo, |argz| < tt - e < vr, F^iz) = Yl'k=i 2 k( 2 k-f)z^>--^ 

Bk, k >1, are Bernoulli numbers. Using this expression it was shown in |Q, 
formula (17), that 


In p{z) 


N , z 
2pq 




z —>■ oo, 


where for arbitrary integer m > 0 we havfl S%{C) = ^m(e) + 0(iV-2—1) 
and the function <hm('C) is defined by the identity 


‘fm(0 


FmiN) - F^(l) 


FmiN-i)-Nr 




®This, of course, follows from estimate of the residual in (122 II see details, e.g. in [S]. 
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withr(r) =Tlni + (l-r)lni^-^(r-p)2, D{t) = ln{l + . 

The function $m(C) is analytic in the union |J^ Bx of unit balls Bx bounded 
by the circles Cx, x E [— j4,j4]. For the function S'^(^) it was shown in |6], 
formula (27), that |5')y(.^)| < for ^ £ Cx and x E [—A, A]. Since we know 
from Theorem [1] that at any point x such that x = Np + [yJ2Npqx) E Z it 
holds that 




M 


E 

u=0 


qujx) I 

Nu/2 I 




therefore it is sufficient to show that 

M 

le—- 




U=0 






for some chosen m = m{M), ^ £ Cx and x £ [—A, A] and some polynomials 
{Pu{O}i^0- ii' i® already shown, then making N large enough we see 
that PuiO — ^ 12 < M, for all ^ due to their coincidence in at 

least [y/N] points. Analyzing functions Fm{i),r{^), D{-^) as functions of 
variable and parameter x we see that for each of these functions we 

I I CX) 

can write the Taylor serieq^ as —)• 0 of the form: ^ for 

j=0 

some integer h. Since x £ [—A, A] we can truncate each series and write 


m 

E uniformly for ^ E {JxBx- It 

is sufficient to set m to be equal to 2M. In the same manner truncating 

m _ , . 

the Taylor series for the function we see that ^ = 

o( jV-M /2 )i for Po(0 = 1 and some polynomials Pu{^), where ^ £ Cx and 
X £ [—A, A], 

Since (|) equals to S^{^) — $m(0 + o{N~'^) we get that \S ^(^)| < 
o{N~^P^) and using the Lagrange Theorem we obtain that — l| = 

o{N-^/^), for C£Cxandx£ [-A, A], □ 


Remark 2. It is possible to obtain estimate le* — X] ^i /2 I = o{N without 

usage of Theorem 1 but directly from the analysis of the function S% at least for small 
values of M. For example, we can obtain the following asymptotic expansions, N —^ oo: 


r 



y/2 (2p-l) X® 

3 y(rT^iV3/2 


+ 0{N-^), 



^/2 {2p - 1) X 


+ 0{N-^). 


^'^We present initial terms of these series in the Remark 2 below. 
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For the function Fm with m = 2 we have: F 2 {z) = 1/ (122:). Analogously we obtain 
expansions: 

TSf? + ° ('"■“'’) ’ + O (n-'=) , 

FiiN) = 1/12 N~^ + 

Therefore we can write the first term of the asymptotic expansion of the function 
taking only the initial terms of the above expressions: 

<i.mW ^ (l-2p) 12x p/ 1 \ ^ qijx) . p/ 1 \ 

23/2(pq)i/2 31^ \N/ \nJ' 

This gives an explicit proof of Lemma 1 for M = 1. It is an interesting question whether 
there is any short and explicit proof of Lemma 1 for all values of parameter M that does 
not use Theorem 1? 

Lemma 2. Let M be a nonnegative integer and A be a positive real. Then 
for X G [—A, A\ and N large enough we have the uniform estimate: 

\df\=0{l). 

a 



Proof. For a given value of M, all x G [—yl,yl] and N = N{A,M) large 
enough we have | I ^ 2 . 


v=0 


2tt 


(x — 0^+^ 


\df\ <2 


0 

and constant C does not depend on x. □ 

Combining Lemma 1 and Lemma 2 we can estimate the second integral 


s! 

_ j 

r - l)df 

27ri J 

(x — 






for ^ G Cx, X G [— 21 , 2 I], that finishes the proof. 
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